基于特征正交分解的一类瞬态非线性热传导问题的新型快速分析方法

提出了一种基于特征正交分解(POD)和有限元法的瞬态非线性热传导问题的模型降阶快速分析方法, 建立了导热系数随温度变化的一类瞬态非线性热传导问题有限元格式的POD降阶模型. 在隐式时间推进方法的基础上有效结合单元预转换方法和多级线性化方法发展了一种加速求解瞬态非线性热传导降阶模型的新型计算方法,并通过二维和三维算例验证了该方法的准确性和高效性. 研究结果表明: (1)降阶模型解的均方根误差在经过初始时段轻微的脉动后稳定于0.01%以下, 而其计算效率比有限元全阶模型提高2$\sim $3个数量级, 并且自由度数量(DOFs)愈大提高的幅度也愈加显著; (2)新型算法解决了常规算法在计算非线性降阶模型时加速性能差的问题, 即使是在DOFs比较小的时候也能够明显提高计算效率; (3)常数边界条件下得到的POD模态可以用来建立相同求解域在各种复杂时变边界条件下的瞬态非线性热传导降阶模型, 并对其传热过程和温度场进行快速准确的分析与预测, 具有很好的工程应用价值.      

A hybrid numerical approach to predict the vibrational responses of panels excited by a turbulent boundary layer

In this work, a hybrid numerical approach to predict the vibrational responses of planar structures excited by a turbulent boundary layer is presented. The approach combines an uncorrelated wall plane wave technique with the finite element method. The wall pressure field induced by a turbulent boundary layer is obtained as a set of uncorrelated wall pressure plane waves. The amplitude of these plane waves are determined from the cross spectrum density function of the wall pressure field given either by empirical models from literature or from experimental data. The response of the planar structure subject to a turbulent boundary layer excitation is then obtained from an ensemble average of the different realizations. The numerical technique is computationally efficient as it rapidly converges using a small number of realizations. To demonstrate the method, the vibrational responses of two panels with simply supported or clamped boundary conditions and excited by a turbulent flow are considered. In the case study comprising a plate with simply supported boundary conditions, an analytical solution is employed for verification of the method. For both cases studies, numerical results from the hybrid approach are compared with experimental data measured in two different anechoic wind tunnels.     

受限亚音速气流中倒置悬臂壁板静气弹稳定性的理论及实验研究

板壳结构在航空航天、高速列车、能量采集等诸多工程领域已经得到了广泛应用. 将悬臂壁板倒置于轴向气流中并在壁板周围流场中设置刚性壁面可有效地调控壁板的失稳速度, 是俘能器优化设计的重要措施之一. 但针对刚性壁面作用下亚音速气流中倒置悬臂壁板的失稳机制仍需要开展深入研究. 本文以受限亚音速气流中倒置的二维悬臂壁板为对象, 以理论分析及风洞实验为手段, 研究了单侧刚性壁面效应对倒置悬臂壁板静态失稳特性的影响规律. 在理论分析中, 首先应用镜像函数法来处理壁面约束条件, 基于算子理论研究获得了以Possio积分方程为表征的壁板气动力, 壁面效应实际表征为一包含移位Tricomi算子的复合算子; 然后将壁板失稳方程的求解问题转化为定区间上的函数逼近问题; 最后, 依据Wererstrass定理并利用最小二乘法求解该最优函数, 以获得系统的失稳临界参数. 在试验研究中依据压杆稳定原理设计了壁板静态失稳的测试方法并完成了风洞实验. 理论分析结果表明, 壁板会发生发散(静气动弹性)失稳, 临界动压随壁板与壁面间距的增加而增大并最终趋于稳定(无壁面情况); 通过理论与风洞实验结果的对比分析, 验证了本文气动力及理论分析的适用性及准确性. 针对倒置悬臂壁板结构的气动弹性失稳问题, 本文提出的方法不涉及系统方程的离散及特征值求解问题, 而是将其转化为了定区间上的函数逼近问题进行求解, 这为弹性结构静气动弹性失稳问题的研究提供了一个可行的新思路.      

Residual stress in polyethylene pipes

This paper deals with residual stress in polyethylene (PE) pipes as an important factor that influences their lifetime. Residual hoop stress distribution in a PE pipe was determined using a methodology previously carried out by the authors on polypropylene pipes. Axial residual stress magnitude was determined by comparison of experimental data and finite element modelling. Based on the obtained results, a new simplified methodology for determining the residual hoop stress is proposed. The method needs only one circular specimen made of pipe, but, unlike older methods, provides a more precise estimation of residual hoop stress distribution taking into account its exponential shape. Some older results from literature are recalculated using this method and residual hoop stress of various PE pipe grades and dimensions are then compared. To provide an idea of residual hoop stress influence on lifetime of a pipe, a lifetime estimation is carried out for the examined pipe.     

Lattice Boltzmann simulation of liquid–vapor system by incorporating a surface tension term

In this study,we investigate the pseudopotential multiphase model of lattice Boltzmann method(LBM)and incorporate a surface tension term to implement the particle interaction force.By using the Carnahan–Starling(CS)equation of state(EOS)with a proper critical pressure–density ratio,a density ratio over 160000 is obtained with satisfactory numerical stability.The added surface tension term offers a flexible choice to adjust the surface tension strength.Numerical tests of the Laplace rule are conducted,proving that smaller spurious velocity and better numerical stability can be acquired as the surface tension becomes stronger.Moreover,by wall adhesion and heterogeneous cavitation tests,the surface tension term shows its practical application in dealing with problems in which the surface tension plays an important role.     

Development and validation of a bioanalytical method by LC‐MS/MS for the quantification of the LAFIS 10 – an antimalarial candidate – and its pharmacokinetics first evaluation

A rapid and highly sensitive method by LC‐MS/MS was developed and validated for the quantification of an antimalarial candidate (LAFIS10) in rat plasma using dexamethasone as internal standard (IS). The chromatographic separation was performed with a Poroshell 120 EC‐C₁₈ column. The mobile phase consisted of water (A) and acetonitrile (B), both containing 10 m m of ammonium formate and 0.1% formic acid, delivered in the form of elution gradient. The LAFIS10 was monitored using an electrospray ionization interface operating in the positive mode in multiple reaction monitoring mode, monitoring the transitions 681.47 → 538.2 for LAFIS10 and 393.20 → 355.30 for the IS. The flow rate was 500 μL/min. The column temperature was kept at 40 °C and the injection volume was 2 μL. The lower limit of quantification was of 10 ng/mL and linearity between 10 and 1000 ng/mL was observed, with an R² > 0.99. The accuracy of the method was >90%. The relative standard deviations intra‐ and interday were <8.80 and <6.37%, respectively. The method showed sensitivity, linearity, precision, accuracy and selectivity required to quantify LAFIS 10 in preclinical pharmacokinetic studies according to criteria established by the US Food and Drug Administration and European Medicines Agency. Copyright © 2014 John Wiley & Sons, Ltd.     

Rapid and sensitive ultra‐high‐pressure liquid chromatography method for quantification of antichagasic benznidazole in plasma: application in a preclinical pharmacokinetic study

Benznidazole (BNZ) and nifurtimox are the only drugs available for treating Chagas disease. In this work, we validated a bioanalytical method for the quantification of BNZ in plasma aimed at improving sensitivity and time of analysis compared with the assays already published. Furthermore, we demonstrated the application of the method in a preclinical pharmacokinetic study after administration of a single oral dose of BNZ in Wistar rats. A Waters® Acquity UHPLC system equipped with a UV–vis detector was employed. The method was established using an Acquity® UHPLC HSS SB C₁₈ protected by an Acquity® UHPLC HSS SB C₁₈ VanGuard guard column and detection at 324 nm_. The mobile phase consisted of ultrapure water–acetonitrile (65:35), and elution was isocratic. The mobile phase flow rate was 0.55 mL/min, the volume of injection was 1 μL, and the run time was just 2 min. The samples were kept at 25°C until injection and the column at 45°C for the chromatographic separation. The sample preparation was performed by a rapid protein precipitation with acetonitrile. The linear concentration range was 0.15–20 µg/mL. The pharmacokinetic parameters of BNZ in rats were determined and the method was considered sensitive, fast and suitable for application in pharmacokinetic studies. Copyright © 2014 John Wiley & Sons, Ltd.     

Modeling and testing of snow penetration

Penetration by a cone into snow is commonly used to characterize snow properties. However, the effects of the diameter and half-angle of the cone on the mechanical properties of snow have not been systematically studied. In addition, no estimation of material parameters in a physically-based model has been made such that the results from penetration provide only an index of snow properties. In this paper, modeling and experimental methods are used to examine the effects of cone geometry on the maximum penetration force and associated hardness, with penetrometers ranging from 2.5 to 4 mm in diameter, 15° to 45° in cone half-angle, and testing both fine-grained and coarse-grained snows. The material point method, in conjunction with the Drucker–Prager cap plasticity model, was used to obtain the theoretical penetration force-distance relationship. Global sensitivity studies were conducted that indicate that the cohesion accounts for 86% of the penetration force, followed distantly by friction angle which accounts for 27%. A general trend, for the simulation results was established: for a given half-angle, the penetration force increases with the increase of diameter which holds for most of the test data as well; for a given diameter, the penetration force decreases with the increase of half-angle, which holds for some of the test data. In addition, for a given half-angle, the hardness decreases with the increase of diameter; for a given diameter, the hardness decreases with the increase of half-angle. To take into consideration the uncertainty of test data, a simple interval-based metric was used to compare test data with simulation results; the comparison was satisfactory. The material parameters from the simulations can thus be considered as calibrated ones for the snow studied.     

膜流动电位测试中Ag-AgCl电极的制备及稳定性研究

膜的动电特性研究中,往往采用流动电位方法。该方法中,Ag-Ag Cl测试电极的质量是影响测试准确性的关键要素。本文采用正交试验设计和方差分析,着重考察电解法制备Ag-Ag Cl电极过程中,电流强度(I)、氯化时间(t)、电解质溶液浓度(CHCl)、电极来源(n)、烘烤温度(T)、活化电极溶液浓度(CKCl)等因素对Ag-Ag Cl电极稳定性的影响。得到Ag-Ag Cl电极的最优制备参数为:电流密度3.0m A/cm2,氯化时间50min,盐酸浓度0.1mol·L-1,烘烤温度120℃,活化电极溶液(KCl)浓度0.001mol·L-1;最显著影响因素为电极来源,显著影响因素为氯化时间和烘烤温度。对最优条件组合下制备的AgAg Cl电极,进行了稳定性实验,结果表明:制备的电极具有较好的稳定性,24h内电极电位漂移量小于0.10mv,7天内电极电位漂移量小于0.2mv;在聚偏氟乙烯中空纤维超滤膜流动电位测试中,流动电位与流动压差具有良好的线性和可重复性,回归曲线R20.99,变化规律符合Helmholtz-Smoluehowski公式,可以较好满足膜的流动电位测试要求。     

A new characteristic approach for incompressible thermo‐flow in Cartesian and non‐Cartesian grids

A virtual‐characteristic approach is developed for thermo‐flow with finite‐volume methodology in which a multidimensional characteristic (MC) scheme is applied along with artificial compressibility. To obtain compatibility equations and pseudo‐characteristics, energy equation is taken into account in the MC scheme. With this inherent upwinding of convective fluxes, no artificial viscosity is required even at high Reynolds numbers. Another remarkable advantage of the MC scheme lies in its faster convergence rate with respect to the averaging scheme that is found to exhibit substantial delays in convergence. As benchmarks, forced and mixed convections in a cavity and in flow over cylinder and between parallel plates are examined for a wide range of Reynolds, Grashof, and Prandtl numbers. The MC and averaging schemes are applied for simulation purposes. Results show the better performance of the MC scheme in forced and mixed convections. Results confirm the robustness of the MC scheme in terms of accuracy and convergence. Copyright © 2015 John Wiley & Sons, Ltd.